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Abstract We investigate the properties of features seen within superluminal sources of- 
ten referred to as components. Our result indicates a fairly strong correlation of r^0.6 for 
quasars, r^0.4 for galaxies, and r^0.8 for BL Lac objects in our sample between compo- 
nent sizes and distances from the stationary core. Assumption of free adiabatic expanding 
plasma enabled us to constrain in general the Lorentz factor for superluminal sources. Our 
estimated Lorentz factor of 7^7 — 17 for quasars, 7^6 — 13 for galaxies and 7^4 — 9 
for BL Lac objects indicate that BL Lac have the lowest range of Lorentz factor. 
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1 INTRODUCTION 

A large number of superluminal motions have been observed in jets of numerous classes of astrophysical 
objects using Very Long Baseline Interferometer (VLB I). The most prominent sources of these jets are 
the Active Galactic Nuclei (AGN), where the jet emission contributes significantly to the spectrum 
of the source (Kording & Falcke 2004). These extragalactic radio jets have been intensively studied 
for a number of important reasons. These include their association with super massive accretion black 
hole, which plays a substantial role in the formation and evolution of galaxies — the building block 
of the universe. Studies and applications of the astrophysics of these radio jets can help reveal black 
holes energetic and their interactions with their environments. Understanding the formation, collimation, 
acceleration, and propagation of these radio jets will aid our understanding of the laws of physics under 
extreme conditions — relativistic speed, high magnetic field strength, etc. 

In the recent years, VLBI observations of AGN jets (e.g. see Kellermann et al. 2004; Britzen et 
al. 2007, 2008; the MOJAVE program; the TANAMI project) have provided us with observations for 
morphological studies and motions of features (often refer to as components) seen within these jets. 
These observations provide us with many values of the apparent luminosity, size, radial distance away 
from the core and proper motion/apparent transverse speed of components moving relativistically along 
the jets. Though these quantities are of considerable interest, the intrinsic properties (Lorentz factor — 
bulk speed/pattern speed, angle to the line of sight) are more fundamental and necessary for kinematic 
characterization of these sources. 

The jets can at the simplest level, be modelled by two intrinsic properties Lorentz factor (7) and 
angle to the line of sight (9). Homan (201 1) noted that by studying the apparent speed distribution of a 
large, flux-density limited sample of AGN jets, we probe the underlying Lorentz factor distribution of 
the parent population, assuming that the apparent speeds of moving jet components are good tracers of 
the underlying flow velocity of the beam (Vermeulen & Cohen 1994; Lister & Marscher 1999). 
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2 THEORY OF RELATIONSHIP 

In the relativistic beam model of extragalactic radio sources (EGRS), the cores of these sources are 
believed to be responsible for the ejection of the observed VLBI jet components, and thus observed 
properties are expected to correlate. As these components are ejected, they move away from the core. 
Assuming ballistic motion, the observed radial distance (D) from the core should be fore-shortened due 
to combined relativistic motion and orientation effect, thus we can write (see Ubachukwu 1998) 

D = D sm9, (1) 

where D is the intrinsic length of the jet component from the core in its rest frame, and 9 the viewing 
angle. The observed apparent speed /3 a with which these components move out from the core from 
orientation based argument, can be defined as 

A, = 7/?<5sin#, (2) 

where 6 is the Doppler factor defined by 6 = 7 _1 (1 — /JcosS) -1 , 7 is the Lorentz factor, which is 
related to /3 the bulk speed defined in units of c (the speed of light) by 7 = ( y/l — /3 2 ) ~ 1 . The maximum 
apparent speed (/3 a ) max for a given Lorentz factor 7 occurs at a critical angle 8 C given by (3 — cos9 c 
which implies that sin6* c = 7 -1 . The Doppler factor for a given Lorentz factor is maximum for (3 = cos# 
and is given by S = 7. Equation (2) can thus be written approximately as, 

(/3a) max «/37- (3) 

Thus, the time (t D ) taken for the component to move away from the core to an observed maximum 
radial distance D max is given by 

*d = ^. (4) 

Kovalev et al. (2005) noted that most new jet features typically increase in size and/or decrease in 
flux density as a result of adiabatic expansion and/or synchrotron losses. Thus, as these components 
move away from the stationary core, they are expected to expand sideways. Assuming free adiabatic 
expansion, we follow (e.g. De Young 2002; Kording & Falcke 2004) to define the expansion speed 
(which is expected to equal the sound speed) as 

ft = , 1 ==, (5) 



+ 
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where m p is the proton mass, n to t total number density of particle, T adiabatic index, P is the ther- 
modynamic equilibrium pressure. The adiabatic index is usually taken to be 4/3 or 5/3 for relativistic 
and non relativistic expansion respectively. Kording & Falcke (2004) showed that for photon emitting 
plasma (which is what we expect from a radio component) the maximal expansion speed will be 

/W = ^r^ 1 (6) 

Thus, the time (i^) taken for the jet component to reach an observed maximum size i? max is given by 

For a free ballistic adiabatic expanding gas in isotropic environment, with constant expansion speed and 
linear speed away from the core, the time to expand to the maximum size £r, should correlate with the 
time to reach the maximum radial distance to - Thus, from equations (4) and (7) we have 
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By implication, a linear regression fit to equation (8) in the form 



logA 



max 



= logi? max + log 




) 



(9) 



will enable us estimate the average Lorentz factor and thus place some useful constraint on the angle to 
the line of sight of the observer. Also, the logD max — logi? max plot is expected to yield a slope of unity, 
which can be tested using a well-defined source sample. We chose a logarithmic form of the relationship 
due to the wide spread in i? max and D max . 

3 DATA ANALYSIS/RESULT 

Our analyses were based on the Caltech Jodrell Bank Flat Spectrum Sources (CJF) defined by Taylor 
et al. (1996) which is a complete flux-limited VLBI sample of 293 flat spectrum radio sources drawn 
from the 6 cm and 20 cm Green Bank Surveys. Of the 293 sources in the original sample, we obtained 
information on 799 jet components for 237 sources from Britzen et al. (2007, 2008), with at least one ob- 
served component apparent velocity for each source, and other parameters like proper motion, distance 
of each component with respect to a stationary core and sizes (major axis) of each component. This 
sample is large, observational strategy and data reduction is homogeneous. Thus, given its size, com- 
pleteness, and range of observed source properties, it will allow a meaningful statistical investigation of 
apparent motions and many correlations of other parameters. 

In our analyses, we assumed that a source contains an ideal beam — one that is straight and narrow, 
with the pattern speed (3 p equal to the bulk speed (3 b . Many sources however, are seen to have more 
than one moving components, with different values of apparent transverse speed (3 a . We selected the 
component with the fastest speed for sources with more than one component, believing that it is the 
component most likely to have the same speed as the beam (see Cohen et al. 2007) and would have 
reached a maximum distance for a given time. More also, Kellermann et al. (2004) suggested that there 
is a characteristic speed for each jet which is probably reflected by the fastest speed observed for the jet 
component. The final sample consists of 177 quasars, 41 galaxies and 19 BL Lac objects. However, 
Britzen et al. (2008) noted that selecting the components based on the brightest components is a good 
representative of the sample and adequate for population studies. Thus, in this paper, we will compare 
our results obtained from selection based on fastest components, with that based on brightest compo- 
nents. We point out that sources with only one observed component, is represented in both selections by 
that component. 

Assuming the radius of the semi-major axis of each source as representing the size, and transforming 
the angular distance of each component from the stationary core to a linear distance, for the selection 
based on fastest/brightest components, we estimated the Lorentz factor of each source from equation 
(8) assuming 4/3 < T < 5/3. The histogram plots of the estimated Lorentz factor for each source in 
each class of object are shown in figures 1 — 6. Generally, the distributions are similar for r = 4/3 and 
r = 5/3 for each class of object, though shifted to higher values for T = 5/3 and are generally positively 
skewed. The median value of Lorentz factor estimated for each class of object for 4/3 < T < 5/3, for 
selection based on fastest components are ^ 9 ± 2 — 11±2;~6±1 — 10 ±2; and ^4±1 — 6±1 
respectively for quasars, galaxies and BL Lac objects. Similarly, for the selection based on brightest 
components we have ^6±1 — 10±2;^5±1 — 8±1 and ^5±1 — 8±1 for quasars, galaxies and 
BL Lac objects respectively for 4/3 < T < 5/3. The quoted errors were estimated from the observed 
error associated with the parameters obtained from Britzen et al. (2007, 2008). 

We also fitted the observed D— R data to equation (9). The result gives : log D = (0.7±0.3)logi?+ 
1.2±0.1 for quasars; \ogD = (0.3±0.5)logi? + 1.1±0.1 for galaxies; and logD = (0.6±0.4)logi? + 
0.9±0.1 for BL Lac objects. This corresponds to a distribution in Lorentz factor of: 7 ~ 7 — 17 for 
quasars; 7 <~ 6 — 13 for galaxies; and 7 <~ 4 — 9 for BL Lac objects for |<T<|, with correlation 
coefficient between D — R as r ~ 0.6, 0.4, and 0.8 respectively; a fairly strong correlation, especially 
for quasars and BL Lac objects. The spread in the Lorentz factor results also imply a distribution in the 
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Fig. 1 Histogram Plot of Estimated Lorentz Factor 
for the selection based on fastest components for 
Quasars (Grey Plot for F = 5/3; Red line for F = 

4/3) 



Fig. 2 Histogram Plot of Estimated Lorentz Factor 
for the selection based on brightest components for 
Quasars (Grey Plot for T = 5/3; Red line for T = 

4/3) 
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Fig. 3 Histogram Plot of Estimated Lorentz Factor 
for the selection based on fastest components for 
Galaxies (Grey Plot for T = 5/3; Red line for F = 

4/3) 



Fig. 4 Histogram Plot of Estimated Lorentz Factor 
for the selection based on brightest components for 
Galaxies (Grey Plot for F = 5/3; Red line for F = 

4/3) 



critical angle to the line of sight 9 C of ~ 3° — 8°, 4° — 10°, and 6° — 16° for quasars, radio galaxies and 
BL Lacs respectively. 

Similar analyses using selection based on brightest components yield: log!? = (0.7±0.3)logi? + 
1.1±0.1 for quasars; logD = (0.9±0.4)logi? + 1.0±0.1 for galaxies; and logD = (0.1±0.4)logi? + 
0.5±0.2 for BL Lac objects, with correlation coefficients r ~ 0.6, 0.5 and 0.3 respectively. The cor- 
responding Lorentz factor for different classes of sources are: 7 ~ 5—12 for quasars; 7^5 — 9 
for galaxies; and 7^2 — 5 for BL Lac objects for |<T<|, corresponding to a distribution in 6 C of 
~ 5° — 10°;^ 6° — 12°; and ~ 12° — 36° for quasars, radio galaxies and BL Lac objects respec- 
tively. The plot of logZ3 max against logi? max for the fastest components is shown in figure 7, while the 
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Fig. 5 Histogram Plot of Estimated Lorentz Factor 
for the selection based on fastest components for 
BL Lac (Grey Plot for r = 5/3; Red line for V = 

4/3) 



Fig. 6 Histogram Plot of Estimated Lorentz Factor 
assuming for the selection based on brightest com- 
ponents for BL Lac (Grey Plot for T = 5/3; Red 
line for T = 4/3) 



plot of log_D max against logi? max for the brightest components is shown in figure 8. In general, the 
Lorentz factor obtained using the fastest components is greater than that obtained using the brightest 
components. We emphasis that values estimated using the selection based on brightest components may 
represent the lower limits to the Lorentz factors for each class of object, and consequently the upper 
limit to the viewing angle. 



4 DISCUSSION AND CONCLUSION 

We have used the relationship between the observed component radial distance D and component size 
R to estimate the Lorentz factor (7) and the corresponding angle to the line of sight 6 for a selection 
based on fastest/brightest components. Our results show a stronger D — R correlation for selection 
based on fastest components than for that based on brightest components for quasars, radio galaxies 
and BL Lac sub-samples. Our results also indicate lack of apparent correlation between D and R in 
the BL Lac data for selection based on the brightest component. This suggests that selection based 
on the brightest components may not be appropriate for the determination of the Lorentz factor or for 
characterizing superluminal motion for this class of object. In general, the result on the analysis of CJF 
sample obtained here, showed that galaxies have lower intrinsic Lorentz factor (lower bulk speed) than 
quasars. There are few BL Lac in the sample, but the result showed that BL Lac have the lowest range 
of Lorentz factor in agreement with (e.g. Gabuzda et al. 1994; Morganti et al. 1995; Kellermann et al. 
2007). This is an indication that BL Lac objects are low power sources (Bicknell 1994). We note that 
equation (9) yields a slope of unity for the plot of the R — D data. The results from our regression 
analyses approximate the theoretical value, especially for quasars, and in general for selection based on 
fastest components. This is an indication that the underlying assumption is plausible. These results do 
show that selections based on fastest component are a good indicator of the bulk speed. 

We note that our estimated values for angle to the line of sight especially for sample selected based 
on brightest components (especially for BL Lacs) are high (see Hovatta et al. 2009; Yang et al. 2012), 
however, Gopal-Krishna et al. (2006) and Wiita et al. (2008) pointed out that allowing for the conical 
shape of ultra relativistic blazar jets with opening angles of a few degrees on parsec-scales, they showed 
that the actual viewing angles of these conical jets from the line-of-sight can be substantially larger 
than the values usually inferred by combining their flux-variability and proper motion measurements 
(Jorstad et al. 2005). This jet geometry also implies that de-projected jet opening angles will typically 
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Log R 

Fig. 7 Plot of logD max against logR max for the fastest components with error bar in standard 
deviation (Legend:Quasars - red diamond; Galaxies - blue circle; BL Lacs - green square) 

be significantly underestimated from VLBI measurements. Also, our estimated Lorentz factor in general 
are lower than that shown in Table 3 of Hovatta et al. (2009), except for galaxies, though their sample 
was obtained at 22 and 37 GHz as against the CJF sample observed at 5 GHz. Kellermann et al. 
(2004) noted that there is a systematic decrease in apparent velocity with increasing wavelength which 
indicates that observation at different wavelengths may be sampling different part of the jet structure. 
Moreover, observations at lower frequency have lower angular resolution and thus more sensitive to the 
lower surface brightness structure located downstream. 

Further analysis of our result, showed that for the sample selected based on brightest components, 
the apparent speed, component size, and component distance from the assumed stationary core are on 
the average less than that of sample selected based on fastest components, but with brighter apparent 
luminosity (see Table 1). They may be expected to have smaller and apparently brighter components, 
but if their motion is unimpeded, they should also on the average, have higher apparent speed. This is an 
indication that their apparent luminosity might be more of projection effect/motion through bends than 
being relativistically Doppler boosted. Thus, Lorentz factor calculated based on brightest components 
for a given sample, will most likely represent the lower limit to the bulk expansion speed, and upper 
limit to the viewing angle, while that estimated from the fastest components will represent the most 
probable values (Gopal-Krishna et al. 2006; Marshall et al. 201 1). In Hovatta et al. (2009), their Figure 
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Fig. 8 Plot of logD max against logR max for the brightest components with error bar in stan- 
dard deviation (Legend:Quasars - red diamond; Galaxies - blue circle; BL Lacs - green square) 

11 shows that jsm8 distribution for their sources peaks around 1, though a number of sources have 
7sin# larger or smaller, thus the assumption of sinff = 1/7 seems reasonable. 

Table 1 Comparison of the Mean values of some parameters. (Between sample selected based 
on fastest component and that selected based on brightest components) 







A, 


Component size 

(mas) 


Component Luminosity 


Component Distance 
(pc) 


Quasars 


fastest 


9.43 


1.14 


4.16 


24.59 




brightest 


5.40 


0.69 


6.87 


13.65 


Galaxies 


fastest 


4.43 


1.30 


1.45 


25.26 




brightest 


2.63 


0.95 


2.28 


20.00 


BL Lacs 


fastest 


4.52 


1.35 


2.40 


10.09 




brightest 


3.45 


0.44 


3.33 


3.73 



Cohen et al. (1988) used a sample of compact radio quasars with measured redshift (z) and proper 
motion /.i to obtained a result which indicated that assumption of relativistic beaming model gives a 
reasonable fit to their data with an upper bound of bulk Lorentz factor lying in the range 9 < 7 < 18. 
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Pelletier & Roland (1989) believed superluminal motion in powerful radio sources can be explained 
by a two fluid model. In their modelling they estimated that relativistic electrons responsible for the 
synchrotron emission of VLBI jets are moving with Lorentz factor 3 < 7 < 10. In the analysis of flux 
variations in BL Lac objects, Mutel (1992) indicated that the Lorentz factor lies in the range 2 < 7 < 4. 
For a sample of 100 sources with published VLBI measurements of the core angular size, Ghisellini 
et al. (1993) assuming Synchrotron Self Compton (SSC) formalism estimated the Doppler factor of 
these sources independent of superluminal motion. Using the estimated Doppler factor and observed 
superluminal speed /3 a they constrained 7 = 10. In VLBI observations of 12 Flat Spectrum Radio 
Quasars (FSRQ) in 2 Jy sample, Vermeulen & Cohen (1994) obtained a distribution of 7 for the fastest 
superluminal components in FSRQ in the range of 5 — 35. Urry & Padovani (1995) showed that for 
their derived luminosity function of FSRQ to match the observed FSRQ luminosity function requires a 
distribution in Lorentz factor in the range of 5 < 7 < 40 with a mean value of (7) = 11. While for 
BL Lac objects, a range of 2 < 7 < 20 and an average of (7) = 3 was needed for their luminosity 
function to match the observed value. They also obtained a critical angle value for maximum beaming 
of 9 C — 14° for FSRQ and 6 C = 12° or 9 C = 19° for BL Lac depending on the assumed model. 
These results generally, agree with our result, especially the result obtained with sample selected based 
on fastest components. The seemingly contradictory result for the BL Lac sample based on brightest 
component may be interpreted as motion through bends with increase in brightness due to projection 
effects. The data are however too poor for any definite conclusion. 

Miller-Jones et al. (2004) pointed out that for tangled magnetic field BocR^ 1 , thus with both syn- 
chrotron loses and adiabatic loses due to expansion, as the component moves out, size increases, thus 
we expect component luminosity to anti-correlate with both component size and component distance, 
except modified by projection/environmental effects. Homan et al. (2002) found that flux density of jet 
components fades from the core as D~ 13 . These are in support of our argument that time for compo- 
nent expansion should correlate with period the components have move a radial distance away from the 
core. The correlation analysis between component size R and component luminosity L shows R — L~ x , 
for the selection based on fastest components with correlation coefficient results are r = —0.3, —0.2 
and —0.3 for quasars, galaxies and BL Lac respectively. These results indicates that selections based on 
fastest component are good indicator of the bulk speed, though without any high level of significance. 

We have shown that expected increase in size as jet components move away from the core can be 
used to characterize the kinematics of these sources. This is true for free expanding jets beyond the 
collimation point. Jets observed at high frequency may be sampling the upstream part (see Kellermann 
et al. 2004), thus a re-analysis of VLBI sources observed at different frequency range may be used to 
place a limit on the point of collimation. 
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